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A LIPSCHITZ METRIC FOR CONSERVATIVE SOLUTIONS OF 
THE TWO-COMPONENT HUNTER SAXTON SYSTEM 

ANDERS NORDLI 


Abstract. We establish the existence of conservative solutions of the ini¬ 
tial value problem of the two-component Hunter—Saxton system on the line. 
Furthermore we investigate the stability of these solutions by constructing a 
Lipschitz metric. 


1. Introduction 


The two-component Hunter-Saxton system, given by 


(l.la) 

(1.1b) 


Ut{x, t) + u{x, t)u^{x, i) = 1 ^ y {u^{z, tY + p{z, <)^) dz 

- j {u^{z,tf + p{z,tf)d^, 
Pt{x,t) + {u{x,t)p{x,t))^ = 0, 


was derived by Pavlov as a model of non-dissipative dark matter |10j . It can also 
be viewed as a high frequency limit of the two-component Camassa-Holm system 
describing water waves [12]. The system (11.11) is a generalization of the Hunter- 
Saxton equation 

( 1 - 2 ) (ut + uUoc)x = 

introduced by Hunter and Saxton as a model of the director field of a nematic liquid 
crystal [5]. 

Here we prove global existence of conservative weak solutions of the initial value 
problem for (ED on the line, and construct a metric that renders the flow Lipschitz 
continuous. Previously Wunsch has proven existence of solutions of ED in the 
periodic setting HD HD, and global existence of dissipative solutions on the real 
line m- Dissipative and conservative solutions are two distinct ways to extend 
the solution past the time where classical solutions break down. Before we discuss 
the difference between conservative and dissipative solutions we will look at how 
solutions break down. 

A common feature for weak solutions of ED and the Hunter-Saxton equation 
[8] is that weak solutions may experience wave breaking, which means that Ux 
tends pointwise to — oo in finite time while u stays continuous. The phenomenon 
is illustrated in the following example. 
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u{x,t) 



Figure 1. A plot of u in Example fOI for t = 0.1.1.5.1.9. 2. 


Example 1.1. Let t S [0,2) and let the functions u and p he defined by 


u(x, t) 


-^t + 1, x<-^t^ + t-l, 

+ t- 1 < X < 0, 
j^x, 0 < X < ^t^ + 1, 


p(x,t) 



X < 0, 

0 < X < ^t^ + 1, 
jt^ + 1 < X. 


Then {u,p) is a weak solution of p.l|) for t G [0,2). Note that Ux{0,t) —>■ —oo as 
t —>■ 2“, which in particular means that wave breaking occurs. We can define the 
energy of the system at time t to be given by 


(1.3) 


f (ulixfi) + p^ix,t)) dx = 2, 
Jr 


which is constant in time, even up to the point t = 2. The energy contained in the 
interval—jt^+t—l < x < 0, given by P^)dx = 1, is also conserved. 

Thus a finite amount of energy is being concentrated in a single point as t —>■ 2~. 


As seen in Example II.II a part of the energy + p^) dx is focused at a single 

point at wave breaking. This illustrates that the energy density is not absolutely 
continuous, but a finite Radon measure in general. Nevertheless, the total energy 
remains constant in time as t —>■ 2~ . Hence Ux,p stay in L^(]R) even if tends 
to minus infinity. This means that the energy can be described by the cumulative 
distribution function of a finite Radon measure. One can extend local solutions to 
global solutions by manipulating the concentrated energy at wave breaking. There 
are at least two ways to extend the solution to a global one past wave breaking. 
On the one hand one could ignore the part of the energy that concentrates on a 
set of measure zero in the continuation, which yields dissipative solutions. On the 
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other hand one could continue by letting the concentrated energy back into the 
system, which would give conservative solutions. In practice that would amount to 
defining u and p by the formulas in Example 11.11 even for t > 2. Thus it is essential 
to include the energy variable in our sets of variables, when constructing global 
conservative solutions. 

We are going to solve the system by the generalized method of character¬ 
istics. The approach is similar to the one by Dafermos in [3], where uniqueness of 
dissipative weak solutions of the Hunter-Saxton equation has been established. To 
that end we map our Eulerian coordinates {u, p, (zi^-l-p^)da;) to the Lagrangian vari¬ 
ables {y,U,H,r) defined as follows. Let y{^,t) be defined by yt{C,t) = u{y{^,t),t), 

= u[y{^,t),t), and define {ux{x,t)'^ + p{x,tY) dx as the 

energy to the left of y{^,t). Note that -I- p'^)dx could be a singular measure. 
We introduce = p{y{^,t),t)y^{^,t). In Section [2] we will rigorously define 

the proper space for the variables (y, U, H, r), and establish mappings between that 
space and the space for conservative solutions of dm. In the above variables the 
Hunter-Saxton system reduces to 


(1.4a) 

yt = u, 

(1.4b) 


(1.4c) 

Ht=0, 

(1.4d) 

n = 0, 

where iJoo 
tion law 

= lim^_ioo iJ(^,0). The time evolution of H follows from the conserva- 

(1.5) 

{ul + p% + {uiul + p^))x = 0, 


see for instance [iin] for the similar conservation law for the Hunter-Saxton equa¬ 
tion. In Section [3] we solve dUD, and together with the mappings from Section [3] 
prove that the we can construct global conservative solutions of dnD. 

In Section[4]we construct a Lipschitz metric. The idea is to construct the metric 
in the transformed variables {y, U, H, r) to avoid having to deal with the measure. 
The Eulerian variables are one-to-one to equivalence classes of Lagrangian variables. 
We construct a functional J that respects the equivalence structure such that it can 
be used as a building block of the metric. The approach here is thus more similar 
to the one employed for the Camassa-Holm equation [316] than the methods pre¬ 
viously used to construct Lipschitz metrics for the scalar Hunter-Saxton equation 

mm- 


2. Mappings between Eulerian and Lagrangian coordinates 

In this section we define the sets of Lagrangian and Eulerian coordinates, and 
investigate the mappings between them. We introduce first an important ambient 
vector space B. 

Definition 2.1. Let Ei be the vector space defined by 

(2.1) Ei = {f G L°°{m.)\ f G L^{R) and lim f{x) = 0}, 

X—>• —OO 

equipped with the norm ||/|| El = ll/lloo + II/'II 2, and E 2 be defined by 

(2.2) E 2 = {fG L°°(R) I /' G L ^ R )}, 
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equipped with the norm 11 / 11^2 = ll/lloo + II/'II 2 - Then define the normed space B 
by B = E 2 X E 2 X El X with the norm 

( 2 - 3 ) ll(/i, /2, h, / 4 )||b = II/1IU2 + II/2IIB2 + II/3II-E1 + Wfih- 

The natural space to look for solutions in Eulerian variables is the following. 
Definition 2.2. The space B consists of all triples (u, p, fi) such that 

(i) UGE 2 , 

(ii) pGL^(IR), 

(Hi) 

(iv) Pac = {ul + p^) dx, 

where Af^(K) denotes the set of positive, finite Radon measures on M. 

We are now ready to define the Lagrangian coordinates as a subset of B. The 
definition is similar to [21 Definition 2.2]. 

Definition 2.3. The set T consists of all quadruples X = {y,U,H,r) such that 
(y — id, U,H,r) G B, and there exists a number c > 0 such that 

(i) y-id,U,H eW^’°°{R),r £ 

(ii) y^ > 0, > 0, + y^ > c > 0 a.e., 

(Hi) y^H^ = f7| + a.e. 

We define the subset To of JF by 

(2.4) Eo = {X£E\y + H = id}. 

We will use the notation Hoo = lim^_>.oo(0 = ll-^lloo- To be able to work 
with the space E we need a mapping from Eulerian to Lagrangian variables, and 
vice versa. Before the mappings are introduced, we will state a useful lemma on 
monotone Lipschitz continuous functions. 

Lemma 2.4 (|7l Lemma 3.9]). Let / : R —>■ R 6e an increasing Lipschitz continuous 
function. Then for any set B with m(i?) = 0, we have that /{ = 0 almost everywhere 
in f~^{B). 

The map L in the following definition maps B into E, and thus represents a 
way to pass from Eulerian to Lagrangian coordinates in a rigorous manner. The 
definition is similar to [U Theorem 4.9]. 

Definition 2.5. Let the mapping L : B ^ Eq be defined by L(u, p, p) = {y, U, H, r) 


where 


(2.5a) 

y{f) = sup{x|^((-oo, x)) + X < ^}, 

(2.5b) 

1 

II 

(2.5c) 

0 

II 

(2.5d) 

»'(0 = {p°y{Ci)ydO- 


Proposition 2.6. The mapping L from Definition \2.5\ is well defined. 
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Proof. The proof follows closely those of [71 Theorem 3.8] and [U Theorem 4.9], but 
here the spaces V and T are different. The difference is that in our case we have 
u,U G L°°(]R), while jU Theorem 4.9] uses u,U G L^(K). To prove that L is well 
defined, let {u, p, fi) G V and define X = {y,U,H,r) = L{u,p,p). We only prove 
that U is in Since u G and y and U are well defined functions, it 

holds that U = uoy is bounded by ||m||oo- D 

We compute an example to illustrate how the mapping L works. 

Example 2.7. Let {u,p,p) G V be defined by 

fl, x<-l, 

(2.6) u{x) = < —X, —1 < a; < 0, 

[o, 0 < X, 

(2-7) p(x) = l[o.i](x), 

(2.8) fJ. = Pac + -^Sq. 

The distribution function F of the measure p is given by 


(2.9) 


F{x) = p{{-oo,x]) 


0, 

X < —1, 

X + 1, 

— 1 < X < 0, 

X + 2 ; 

0 < X < 1, 

5 

2> 

1 < X. 


Then {y, U, FI, r) = L{u, p, p) is given by 


(2.10a) 


(2.10b) 


(2.10c) 


(2.10d) 


yii) 


'C, ^<-1, 

|(e-i), 

< 0, 

i(£ —3) 3 ^ ^ 7 

2V> 2 2’ 2—^ — 2’ 

2 ’ 2 — 


u{0 


H{f) 




[l, e<-l 

^-i(e-l), -1<?<1, 

[o, i<e, 

'o, e<-l, 

Ke + I), 

< e, 1 <^< I, 


See Figure\^ for a plot of the functions F, y, and H. Note that the Dirac delta 
corresponds to a flat interval in y, with the length of the interval equal to the strength 
of the delta. 


We also need a mapping M : F ^ V that takes us back from Lagrangian to 
Eulerian variables. The definition is similar to [H Theorem 4.10]. 
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Figure 2. A plot of the functions F, y, and H in Example 12.71 
Note that the jump in F corresponds to the flat interval in y. 

Definition 2.8. Let the mapping M : F ^ V be defined by M{y, U, F[, r) = {u, p, p) 
where 

(2.11a) u{x) = U{y{^)), for some ^ such that x = J/(C)) 

(2.11b) pda: = y#(rdO, 

(2.11c) p = y^(iJjdO- 

The notation /#(/r) denotes the push-forward of the measure p, by the measurable 
function f, i.e. f^(p){A) = dp for all measurable sets A. 

Since there are instances where y^ = 0 on some interval, see Example 12.71 one 
might encounter difficulties when trying to invert y. Also it is not clear that the 
range of M is F. It is therefore necessary to prove that the mapping is well defined. 

Proposition 2.9. The mapping M in Deftnition \2.8\ is well defined. 

Proof. The proof follows closely those of [3 Theorem 3.11] and [U Theorem 4.10], 
but here the spaces F and F are different. The differences are that here u,U G 
L°°(]R), while [T] Theorem 3.11] and (4] Theorem 4.10] use u,U G L^(]R). Let 
X = {y,U,F[,r) S F and (u, p, p) = M{X). We prove that u G L°“(R) since 
everything else is covered in the proof of [U Theorem 4.10]. Since u is a well 
defined function we must have that u is bounded by ||f7||oo- D 
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Example 2.10. Let X G J-'o be given by L{u,p,g.) in Example \2. 7[ Then M{X) = 
{u,p,p). To see that this is the case we compute the distribution function of the 
measure we get by applying M to X. Let a; G R and M{X) = (it, p, p), then 

p{{-oo,x])= f 

0 , X < — 1 , 

5 d^ = a: + l, —l<a:<0, 

5d^+d + l = a; + |, 0<a;<l, 

- 1 < T 

which is equal to the distribution function F of p from Example \2.7\ plotted in Figure 

H 


In the Lagrangian formulation there are four unknowns, while in the Eulerian 
there are only three. Hence it is not surprising that there is some redundancy 
in the Lagrangian formulation. Indeed there are equivalence classes in J- such 
that all elements in an equivalence class map to the same element in T). The 
equivalence classes are determined by a group of relabeling functions, G, and a 
relabeling operator •. 

Definition 2.11. We define the group G and the group action • of G on E as 
follows. 

(i) Define G as the group of homeomorphisms / : R —> R such that both /—id £ 
ITi’°°(R), /-I - id G ITi’°°(R), and /e - 1 G L^{R). 

(ii) Define a group action • : E x G ^ E by {X, /) i-A {y o f,U o /, iL o /, (r o 
f).f') = X.f. 

The result is that X and X • f \s mapped to the same element in T>. For more 
details on the group G and the group action •, see [H Proposition 4.5] and [3 
Proposition 3.4]. 

Proposition 2.12. Let f G G and X G E, then M{X • /) = XI{X). 

Proof. This proof is similar to that of [3 Theorem 3.11]. Let X G E and / G G be 
given and let {u, p, p) = XI (X), and {u,p,p) = XI{X • f), respectively. For a proof 
that u = u and p = p see [3 Theorem 3.11]. We prove that p da: = p da; in the 
sense of measures. Let H C R be of finite measure, and recall that pda; = yy{rdf). 
Then 


pda;(A) = {y o f)y,{r o f f^ df){A) 

(2.12) = [ ro/(0/c(0de 

d(yo/)-i(A) 

Since / is invertible and Lipschitz continuous, and r G L)^^(R) we can change 
variables to obtain 


/ ro f{^)f^{^)df= ( 

d(yo/)-i(A) Jf 


r(^)dC 


foiyof)~^{A) 

[ r{^)d^ = yy{rd^){A). 

ly-HA) 


(2.13) 
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In the last equality we have used that 

(2.14) /(U I y o f{0 e A}) = {/(O I y o /(C) e A} = {C I y(C) e A}, 

since / is onto, continuous, and strictly increasing. Equation (I2.13|) implies that 
p = p almost everywhere. □ 

We need that when one maps an element of Eulerian coordinates to Lagrangian 
coordinates and back that one should end up with the same element. In lieu of 
the previous propositions the converse cannot hold for F. If we restrict M to 
however, we get that the composition L o M is the identity function on Fq. 

Lemma 2.13 ([3 Theorem 3.12]). The functions L and M satisfy 

(2.15) M o L = id®, 

(2.16) LoM = id^„, 
when M is restricted to Fq. 

3. Existence of solutions 

In the previous section we saw that the space of Eulerian coordinates could be 
represented by Lagrangian coordinates. We now want to reformulate the initial 
value problem of dm in Lagrangian coordinates. In this section we motivate the 
system ca), and show the existence of solutions for this system. Then we define 
conservative weak solutions of (da, and show that we can construct such solutions 
by mapping the initial data from Eulerian to Lagrangian coordinates, solve dm, 
and map the solution back to Eulerian coordinates. 

Due to the fact that a finite amount of energy accumulates in a point at wave 
breaking, we replace (u^ + p^) dx by a measure p such that pac = (w^ + P^) dx. 
Then the conservation law dEU reads 

(3.1) pt + {up)x = 0. 

Let f G R and define y by 

(3-2) ^y{^,t) =u{yif,t),t), yif.,0) =yo{f,). 

Then we define U{^,t) = u{y{^,t),t), and for y{£,,t) G (supp ^s(t))°, we define 
= p{i —oo,y{^,t)]). We assume that y{^,t) G {snpp ps{t))^ and get 

(3.3) ^ Myitt)F) + Ui(„t)u^{y{^,t),t) = ^H{f,t) - 

where C/(^,0) = '«(yo(C))0)j and we have used the equation for u in (11.11) . Since 
y{£,,t) € (supp ^s(t))'^ we have from the conservation law (13.11) that 


dt 


H{f,t) = Pt {{-oo,y{^,t)]) 


dp 

dx 


{y{^,t),t)U{^,t) 


(3.4) 


= -UiU)^iy{^,t),t) + U{t)^{y{^,t),t) = 0, 


where 0) = p[{ —oo,j/o(0])- Since p is a conserved variable it is natural to 
look at p(x) dx = p(2/(0)%(C) Define now r{t) = p{y{t),t)y^{t), then 

d 


dt 


r{t) = iptiy{t),t) + p^{y{t),t)U(t))y^{t) + p{y{t),t)U^(t) 

= -p{y(t), t)uS) + piy(t),t)u^it) = 


(3.5) 
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(3.6) r{0) = p{yo,0)yo^, 

where we have used the equation for p in (HU. The derivation assumed that y(t) 
was outside the support of the singular part of p{t), but we will extend the system 
to all of K.. 

The system ()1.4|) can be solved explicitly. We will be interested in the initial 
value problem with initial values in J^q. 


Proposition 3.1. 

(3.7a) 

(3.7b) 

(3.7c) 

(3.7d) 

in Tq, is given by 


The solution of the system dni) with initial data 

y\t=o = yo, 

U\t=o = Uq, 

H\t=o = Hq, 

r\t=o = "^0, 


(3.8a) y(^,t) = i(77o(0 - + Uo{Ot+ yo{^), 

(3.8b) 17(e,t) = + 

(3.8c) H{C,t)=Ho{0, 

(3.8d) r(^,t)=ro(0, 


where Hoo = ||77||oo = lim^^oo H{f). 

Proof. To find the solutions we integrate (11.41) with respect to t, starting with the 
equations for r and H, and then proceed to U and finally y. Uniqueness follows 
from the linearity of the system. □ 


Example 3.2. Let X = L{u,p,p) be as in Example\2.7\ Then St{X) equals 




U{U) 


~f6^^ + t + ^, 

i(e-i)i^ + i(e-i), 

+ e<-i 

.8^’ 2 — 


-1, 

1<C< f, 

l<e- 






Theorem 3.3. The solution operators St '■ J- ^ J- constitute a semigroup. Fur¬ 
thermore the semigroup is Lipschitz continuous in T in the sense that for X, X G P 
there holds 

||5t(X) - 5 *(X)||b < (it" + t + 1)||X - Xb. 


(3.9) 
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Proof. Let Xq = {yQ,Uo, Ho,ro) denote the initial data and St{Xo) = X{t) = 
the solution of (|1.4I1 at t. We need to show that the solution 
is in X. For each t we have that H is bounded by ||-ffo||oo, \U\ by |||-ffo||oo^ + 
||C/o||oo and \y\ by |||Lfo||ooi^ + ||t/o||ooi+||2/o||oo- Since the solutions (15^ are linear 
combinations of the initial data plus a constant we can differentiate the solutions 
with respect to If we differentiate the solutions (13.81) with respect to ^ we obtain 


(3-lOa) yd^,t) = + Uo^{0t + yod0, 

(3.10b) udi,t) = lHodOt + UodO, 

(3.10c) HdC,t) = HodO- 

We have the following estimates 


r{t) = ro, 

0 < Hdt) < Ho^, 

|blc(0l < ■*" l^osl’ 

(3-11) lUiit) ~ 1| < + \Uo^\t + |?/oc — 1|, 

which are square integrable and bounded. Furthermore, y— id, U,H G 1F^’°°(R) as 
this holds for the initial data, and for each t the solutions are linear combinations 
of the initial data. Thus property (i) in Definition 12.31 is proved. Consider now 

C/| = + Uo^Hod + Uq^, 

2 2 
r =ro, 

(3.12) y^H^ = iiLoV + Uo^Hod + + rl 


where it has been used that Ho^yo^ = + Tq. The above proves that (in) 

in Definition 12.31 holds. Non-negativity of y^ follows from being non¬ 

negative, and thus y^ has to be non-negative due to (13.1211 . We show that there 
exists a c{t) dependent on t such that {y{t) + > c{t) > 0. By assumption it 

holds for t = 0 with a constant c(0), and since for each ^ the functions y^ and 
are continuous in t it will hold on some interval [0,r(^)). We choose T(^) to be 
the maximal time for which it holds. Then for t G [0,T(^)) we have 


d 1 
dt y^ + 


< 


(y? + 

1 \Ui\ 
% + H^y^ + 


< 1 

~ + H^y^ + 

1 1 


(3.13) 
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By Gronwall’s inequality 

(3.14) 


% 




62 % 


where c(0) is the initial constant, and T{^) can be chosen to be arbitrarily big, 
which shows that property {ii) in Definition 12.31 holds. We prove (13.91) . Since the 
time evolution of the derivatives are linear in derivatives of X we can use (13.101) to 
find 

(3.15) Wyd't) ~ yd^)h < - Ho^\\2t'^ + ||Go{ - Uoiht + ||?/0C - yoclb, 

(3.16) ||G^(t) - U^{t)\\2 < - Ho^ht + ||17o{ - Dojlb- 

In the L°“-part the term makes the solution operator nonlinear. However, we 
can bound 

(3.17) \\H-H- + iffoolloo <l\\H- H\\oo, 

and thus 

(3.18) || 2 /(f) - yit)\\oo < ^\\Ho - Ho\\oot‘^ + \\Uo - UoWoot + bo - j/olU, 

(3.19) \\U{t) - i7(t)|U < ||i7o - HoWoot + ||C/o - Holloo. 


Since H and r do not change in t, the estimate is proved. The estimate (13.1711 
implies that the right-hand side of C3D is Lipschitz in X, and thus St satisfies the 
semigroup property. □ 


We define the map Tj : 27 —)• 27 by 
(3.20) Tt = MoStoL. 

Definition 3.4. A triple {u, p, fi) € T) is said to be a conservative weak solution of 
(11.11) if for any test function cj) € C'“(]R x [0, oo)), 

(3.21a) 

OO X OO 

J J [ucft + + J “ J (j) dxdt = - J {ucj)\t=o) da;, 

0 R —OO X R 

(3.21b) 

OO 

J J {pift + pu4>x) dxdt = - J {p(j))\t=o dx, 

OR R 

(3.21c) 

oo 

J J {(j)t U(l)x) dfi{t)dt = J dt, 

R 


0 R 

and in addition 


(3.22) 

holds for all t > 0. 


We can now use the operator T* defined by (13.201) to construct conservative weak 
solutions of dni). 
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Theorem 3.5. The operator Tt : (u, p, fi) i—>■ {u(t), p{t), p(t)) maps an initial value 
to a conservative weak solution of CH) in the sense of Definition\3.4\ 


Proof. The idea of the proof is to first map the initial data in V to J-q by L, then 
solve the problem there and for each t map the solution to V by M. Change 
of variables must be done on the set | y^{^) > 0}, but since both and r 
equals zero almost everywhere on the complement we can integrate over K. when we 
change variables in u and p. Since p(t) = y(t)^ {H^ d^) it holds that p ((—oo, x]) = 
sup{H{^) I y{^) = x}. Thus 

OO a; OO 

J J {u(j)t + 4>X + j dp{t)-J dp{t))(l)dxdt 

0 M —oo X 

oo 

= J J U{(ft o y + o d^dt 

0 E 

oo / OO \ 

+ ^ y y y <iT{t) ~ \ J (l>°yye. d^dt 

0 E y—oo —OO J 

oo 

U{i,t){^(l){y{^,t),t) - ^U(t):^{y{^,t),t)) 


0 E 


- ^ (H{£„t) - ii?oo ] (foy 


yi d^dt 


[^{Ucfoy)- ^u^(j}^ o y]y^ d^dt 


OO 

= y J [^(Uffoyyi) - {^u‘^(foy)i] d^dt 

0 B 

= —y uo{x)(j){x,0) dx. 

R 

The equation for p is treated in the same way, 

OO OO 

y J p(j)t + {pu)(j)x dxdt ^ y y + U(l)x)pdxdt 

{(ft o y + U(fx o y)r d^dt 

0 R 

OO 

y J ■^i4'{y{^,t),t)r{^,t)) d^dt 


0 E 
oo 
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= - J </'(2/o(0)0)ro(0 dC 

= — J 0)poix) dx. 


The equality on conservation of p, is proved 

OO OO 

/ J^ J J+ u(j);^) o yH^d^dt 

0 R 

OO 

= J J H^^{(l)oy) d^dt 


0 R 


0 R 


= - j H^{(j)oy)\t^o dC 

M 

= - J (t)\t=o dp\t=o- 

R 

For every t we have that = R, and thus 

p{t)(R)= J H^{^,t)d^ = J Ho^d^ = pq(R). 


y(t)-i(«) 


□ 


Example 3.6. Let {uo,po,po) € T) be as in Example \2.L\ Then Tt{uo, po, po) is 
given by 


(3.23) 


(3.24) 

(3.25) 



f-ft+l 



X < — +t — 

1, 


1-it 

- H 

8 ’ 

-^t"^ + t- 1 < 

a;<-^t2 

j? 

II 

2x 
t ’ 



_J_f2 < 

16^ — — 16 

2 

5 


1 ^+7 ^ 




+ 1? 


L 8 ’ 



+ 1 < a:, 



[0, 

X < 

1 

16 



p{x,t) = < 

ItVT’ 

16 

< 

1—1 

+ 

VI 

H 



[o, 

5 f2 
16 

+ 

1 < a;, 


p[t) = p{t)ac + ^ 

d(o,o) 

(a; 

,t) + <5(_i^2)(a^:0 



Remark 3.7. Note that in Examvle \l.l\ and \3.6\ the functions u{x, 0) and p(x, 0) 
coincide, while in Example \1.1\ and Example \3.(\ respectively, we have p(0) = 
{ux(x, 0)^ + p{x, 0)^) dx and p{0) = {ux{x, 0)^ + p{x, 0)^) da; + ^(5o. The solutions 
in Example \1.1\ and \3.6\ differ, hence it is important to include the measure p in 
the description of global solutions and the definition of the Lipschitz metric. 
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4. The Lipschitz metric 

In this section we construct a metric on Tq that renders the flow Lipschitz con¬ 
tinuous with respect to initial data. We saw in Proposition 12.121 that the mapping 
M : jF ^ T) is relabeling invariant. An important fact is that relabeling commutes 
with the solution operator St- 

Proposition 4.1. For any X G T, f G G it holds that St{X • f) = St{X) • f. 
Proof. Any component of St{X) is a linear combination of components of A. □ 

The problem of creating a metric directly on X is that X and X may correspond 
to the same solution in Eulerian coordinates even if A A. Thus we will try to 
compare solutions in Aq, but we must prove that we can reach Aq from all of A via 
the group action • from Definition 12.111 (ii). 

Definition 4.2. Define the map II : A —>■ Aq by 

(4.1) nX = X*{y + H)-\ XgT. 

To ease the notation we write IIA, despite the fact that II is not a linear operator. 

For the map II to be a relabeling we need that y + H G G. 

Proposition 4.3. Let X G Aq, and X{t) = St{X). Then for allt >0 we have that 
y{t) + H{t) G G, and e~^* < y^(t) + H^(t) < jt^ -|-t-|- 1 for almost every ^ S K. 

Proof. Let y + H = /, we show that f G G. From the definition of A we have 
that f — id G IT^’°“(R) with f^ — 1 G L^(]R), with c < /j < C for some positive 
numbers c and C. Thus there exists a Lipschitz continuous inverse f~^ such that 
f~^ — id G IT^’°°(R). Hence f G G. The lower bound on y^{t) + H^{t) is given by 
(I3.14|) . while the upper bound is a result of the time evolution and Ij/o^Ij and 

I C/ 05 1 being less than 1 . □ 

The metric induced by || • || b will unfortunately give a positive distance between A 
and A • / , even though they will map to the same element in Eulerian coordinates 
via the mapping M. One could potentially restrict attention to the class Aq by 
comparing X • {y + H)~^ and X • {y + H)~^, but it has proven difficult to control 
the t-dependence of {y + H)~^. Instead we minimize the distance over all possible 
relabelings. Following we define J : A x A —>■ R by 

(4.2) J(A, A) = inf (||A . / - A||b + || A - A . sWb) . 

If we instead tried ||A • f — X • g\\B we would not be able to separate r from —r 
as the next example illustrates. 

Example 4.4. Let r = l[o,i]) and y,U,H such that X = {y,U,H,r) G Aq. Then 
X = {y, U, H, —r) G Aq as well. Let 0 < e < 1 and define f^GGby 

U, f<0, 

no ={e^, 0 < ? < I. 


(4.3) 








LIPSCHITZ METRIC FOR THE HUNTER-SAXTON SYSTEM 


15 


Then for any t>0, all components of St{X) and St{X) except r,r are equal, and 
thus, 


(4.4) 


\\St{X).r 


St{X).f\\B< 



< 2-v/e, 


and hence the infimum over all f,g £ G must equal zero. 


In the way we have defined J here we avoid the scenario in Example 14.41 since 
we cannot make both r and f small at the same time. Now, J will not separate 
X and X • f, and behaves as expected with t. However J is not a metric, as the 
triangle inequality fails. One can salvage a metric by taking the infimum over sums 
in J over finite sequences in 


Definition 4.5. Let X,X ^ Xq, then define d : To x J-q ^ M by 

N 

(4.5) d{X, X) = inf ^ J(X„_i, N„), 

n— 1 

where the infimum is taken over all finite sequences {Xn}n=o such that the 

endpoints Xq and X^ satisfy 

(4.6) X = Xo, 

(4.7) X = Xn. 


It is not at all clear that d only vanishes when X = X. The purpose of the next 
lemma is to assert that we have a positive lower bound on d{X, X) when X differs 
from X. 


Lemma 4.6 (O Lemma 3.2]). For any X,X & To we have 

(4.8) \\y - yjloo + \\U - C7||oo + \\H - i7||oo < 2d(X, X). 

Lemma 14.61 states that if the distance between X and X equals zero, then 
(y, U, H) and {y, U, H) coincide. Still, r and f could, in principle, differ. The 
next lemma shows that this cannot be the case, and consequently d is a metric on 
To. 

Lemma 4.7 (0 A weaker form of Lemma 6.4]). Let X,X be in To, then if 
d{X, X) = 0 we have that r = r. 

We need to estimate j(n5't(A), nS't(A)) in terms of j[St{X), St{X)). 

Lemma 4.8. For X,X G To, there holds 

(4.9) J(n5t(A),n5t(X)) < e^*J{St{X),St{X)). 

Proof. The proof consists of two parts. First we show that for f G G with 
\/ll/^l|oo < C, with G > I, it holds that for any X,X G T 

(4.10) J{X • f,X) <GJ{X,X). 

Second we show that if A G To then {y{t)+L[{t))~^ is in G with {{y{t)+L[{t))~^)^ < 
6 2 *. The L°°(M.) part of J is invariant with respect to relabeling. It suffices to show 


^This idea is due to A. Bressan. 
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that it holds for the L^(K.) part. To that end let h and h belong to T^(K.), and let 
g G G. Then by change of variables and denoted f = go 




/ \hoff^-hogg^\^d( 

JR 

= f \Hv)fdf~\r])) - ho f{g)g^{f-\g))\^ 

JR 

= [ \h{v)-ho f{g)f^{g)ff^o f-^{g)dg 

JR 

(4-11) < ll/dloo / \h{g)-ho f{g)f^{g)\'^dg, 

JR 

where we have used that ■ We get the inequality 

(4.12) inf ||/io//^ -hogg ^\\2 < C inf \\h-hoJJ^\\ 2 , 

g&G f^G 


dg 


with C = ^||/{||oo- The part where one relabels ho ff^ is fine since X* f*g equals 
X • / for some f G G, and infimum over G is taken afterwards. Proposition 14.31 
states that y{t) + H{t) G G, and that < Vi + + t + 1). One can 

then invert the lower bound on and apply the first result twice to obtain 

the bound in the lemma. □ 

We are now ready to prove the Lipschitz theorem on jFq. 

Theorem 4.9. Let X,X G Xq, then for all t > 0 it holds that 


(4.13) 


d(n5t(x),n5,(x)) < e^dld + t + i)d{x,x). 


Proof. Let 1 > ^ > 0 and X,X G be given and choose {Xn}n^Q, {/n}n^i, and 
{9n}d=o such that Xq = X,Xn = X and d{X,X) + e > ||X„ • /„ - X„_i||b + 
||X^ — Xn -1 • gn-i\\B- Then from the definition of d we have 

N 

d(n5t(x),n5t(x)) < ^ j(n^t(x„),n5t(x„_i)) 


N 


(4.14) 


<edY^j[St{X^),St{Xr,-i)), 


n=l 


by Lemma [T?51 From the definition of J we get 


N 


d{USt{X), USt{X)) < e5* ^ J{St{X^), 5t(X„_i)) 


n—1 


N 


(4.15) 


< ei\-r + t+l)J2 ( 11 ^" • fn - 
n—1 

T ll-^n X<fi—i • —ills) 

“ T S'). 


The inequality holds for each s in the range (0,1), which implies that 
(4.16) d (n5,(x(t)), nst{x{t))) < e^d^d +1 + i)dix,x). 
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□ 

Since Tq is in one to one correspondance with T) the metric d on J^q induces a 
metric dx> on P. 

Definition 4.10. Define the metric dx> on V by 

(4.17) d-D{{u,p,fi), iu,p,fl)) = d{L{u,p,p),L{u,p,p.)), 
for any (u, p, p), {u, p, p) S V. 

Theorem 4.11. The solution operator Tt defined by 1)3.201) forms a Lipschitz con¬ 
tinuous semigroup on {D,du) in the sense that for any {u, p, p), {u, p, p) € V the 
inequality 

(4.18) dv{Tt{u, p, p),Tt(u,p,p)) +t+l)d-D{{u,p,p)fiu,p,p)) 

holds. 

Proof. The Lipschitz continuity is a corollary of Theorem [23] and Theorem l4.9l We 
prove the semi group property. From Proposition 12.121 we have M o n = M, hence 

(4.19) Tt = M o Sto L = M oUo Sto L. 

From Theorem 13.31 and Proposition 14.II we have 

TtTs = MonoS'ioLoMoIIoS'soL = M oUoStoUoSsoL = MoIloSt-\-s^L = 

□ 
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Helge Holden, and Xavier Raynaud. 

References 

[1] A. Bressan and A. Constantin. Global solutions of the Hunter—Saxton equation. SIAM J. 
Math. Anal., 37(3):996-1026, 2005. 

[2] A. Bressan, H. Holden, and X. Raynaud. Lipschitz metric for the Hunter—Saxton equation. 
J. Math. Pures AppL, 94(l):68-92, 2010. 

[3] C. M. Dafermos. Generalized characteristics and the Hunter—Saxton equation. J. Hyperbolic 
Differ. Equ., 8(01):159-168, 2011. 

[4] K. Grunert, H. Holden, and X. Raynaud. Global solutions for the two-component Camassa— 
Holm system. Comm. Partial Differential Equations, 37(12):2245-2271, 2012. 

[5] K. Grunert, H. Holden, and X. Raynaud. Lipschitz metric for the Camassa-Holm equation 
on the line. Discrete Contin. Dyn. Syst. Ser. A, 33(7):2809-2827, 2013. 

[6] K. Grunert, H. Holden, and X. Raynaud. Lipschitz metric for the two-component Gamassa- 
Holm system. In Ancona F. et ah, editors, Hyperbolic Problems: Theory, Numerics, Appli¬ 
cations, pages 193-207. American Institute of Mathematical Sciences, 2014. 

[7] H. Holden and X. Raynaud. Global conservative solutions of the Gamassa-Holm equation-A 
Lagrangian point of view. Comm. Partial Differential Equations, 32(10):1511-1549, 2007. 

[8] J. K. Hunter and R. Saxton. Dynamics of director fields. SIAM J. Appl. Math., 51(6):1498— 
1521, 1991. 

[9] J. K. Hunter and Y. Zheng. On a completely integrable nonlinear hyperbolic variational 
equation. Phys. D, 79(2—4):361—386, 1994. 

[10] M. V. Pavlov. The Gurevich—Zybin system. J. Phys. A, 38(17):3823-3841, 2005. 

[11] M. Wunsch. On the Hunter-Saxton system. Discrete Contin. Dyn. Syst. Ser. B, 12(3):647— 
656, 2009. 

[12] M. Wunsch. The generalized Hunter-Saxton system. SIAM J. Math. Anal., 42(3):1286—1304, 

2010 . 






18 


A. NORDLI 


Department of Mathematical Sciences, Norwegian University of Science and Tech¬ 
nology, NO-7491 Trondheim, Norway 
E-mail address: andenord@math.ntnu.no 


